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Abstract. The purpose of this work is to present accurate and 
detailed mass - radius relations for white dwarf (WD) models 
with helium, carbon, oxygen, silicon and iron cores, and with 
and without a hydrogen envelope, by using a fully updated stel- 
lar evolutionary code. We considered masses from 0.15 Mq to 
0.5 Mq for the case of helium core, from 0.45 Mq to 1 .2 Mq 
for carbon, oxygen and silicon cores, and from 0.45 Mq to 1 .0 
Mq for the case of an iron core. In view of the recent measure- 
ments made by Hipparcos that strongly suggest the existence 
of WDs with an iron - dominated core, we focus our attention 
mainly on the finite - temperature, mass - radius relations for 
WD models with iron interiors. In addition, we explore the ef- 
fects of gravitational, chemical and thermal diffusion on low 
mass helium white dwarf models with hydrogen and helium 
envelopes. 

Key words: stars: fundamental parameters - stars: white 
dwarfs - stars: evolution - stars: interiors 



1. Introduction 

It is a well known fact that about 90% of stars will end their 
lives as white dwarf (WD) stars. At present we know differ- 
ent routes that drive stellar objects to such a fate. It is widely 
accepted, for instance, that low mass WDs with stellar masses 
M ^ 0.45 Mq are composed of helium and that they have 
had time enough to evolve to such state as a result of binary 
evolution. For intermediate mass WDs, stellar evolution theory 
predicts an internal composition dominated by carbon and oxy- 
gen. Finally, for the high mass tail of the WD mass distribution, 
theory predicts interiors made up by neon and magnesium. 

Over the years, it has been customary to employ mass - 
radius relations to confront theoretical predictions on the in- 
ternal composition of WDs with observational data. This is so 
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because, as it is well known since Hamada & Salpeter (1961) 
(hereafter HS, see also Shapiro & Teukolsky 1983), zero - tem- 
perature configurations are sensitive to the internal composi- 
tion. One of the effects that allow us to discriminate the WD 
internal composition for a given stellar mass is related to the 
dependence of the non - ideal contributions to the equation of 
state (EOS) of degenerate matter (such as Coulomb interactions 
and Thomas - Fermi corrections) on the chemical composition. 
These contributions to the EOS are larger the higher the atomic 
number Z of the chemical constituent. Another very important 
effect is that, in the case of heavy elements like iron, nuclei 
are no longer symmetric (Z = 26, A = 56 for iron), yielding 
a mean molecular weight per electron higher than 2. Accord- 
ingly, for a fixed mass value, the WD radius is a decreasing 
function of Z. 

Recently, Provencal et al. ( 1 998) (and other references cited 
therein) have presented the Hipparcos parallaxes for a hand- 
ful of WDs. These parallaxes have enabled to significantly im- 
prove the mass and radius determination of some WDs, thus 
allowing for a direct confrontation with the predictions of WD 
theory. In particular, the suspicion that some WDs would fall 
on the zero - temperature, mass - radius relation consistent with 
iron cores (see Koester & Chanmugam 1990) has been placed 
on a firm observational ground by these satellite - based mea- 
surements (see Provencal et al. 1998) (however, see below). In- 
deed, some WDs have much smaller radii than expected if their 
interior were made of carbon and oxygen, suggesting that, at 
least, two of the observed WDs have iron - rich cores. Specif- 
ically, the present determinations indicate that Procyon B and 
EG 50 have radii and masses consistent with zero - tempera- 
ture, iron WDs. Obviously, such results are in strong contradic- 
tion with the standard predictions of stellar evolutionary calcu- 
lations, which allow for an iron - rich interior only in the case 
of presupernova objects. Although these conclusions are based 
on the HS zero - temperature, mass - radius relations (note that 
EG 50 has an effective temperature. Toff, of Tcff^ 21000 K), it 
is clear that unless observational determinations are incorrect, 
the interior of the above - mentioned WDs is much denser than 
expected before. 

Before the above - mentioned determinations, an iron com- 
position has been considered as quite unexpected. In fact, the 
only attempt of proposing a physical process able to account 
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for the formation of iron WDs is, to our knowledge, that of Is- 
ern et al. (1991). In their calculations, Isern et al. find that an 
explosive ignition of electron - degenerate ONeMg cores may, 
depending critically upon the ignition density and the veloc- 
ity of the burning front, give rise to the formation of neutron 
stars, thermonuclear supernovae or iron WDs. It is therefore 
not surprising that, apart from the study carried out long ago 
by Savedoff et al. (1969), who did not consider the effects of 
electrostatic corrections, convection and crystallization in their 
calculations, very little attention has been paid to the study of 
the evolution of iron WDs. 

We should warn the reader that the existence of WDs with 
an iron - rich interior is still under debate. In particular, despite 
recent claims of an iron - rich interior for Procyon B, in the 
report of this work our referee has told us about a new reanalisis 
of the observational data which, in a preliminar stage, seems to 
indicate an interior composition for this object consistent with 
a carbon one. 

Another interesting possibility is that these objects may 
content some extremely compact core, as proposed by Glen- 
denning et al. (1995a, b). They suggested the existence of stel- 
lar objects composed by a strange quark matter (with a density 
of « 5 X 10^^ g cm^'^) suiTounded by an extended, normal mat- 
ter envelope. These configurations have been called "strange 
dwarfs". It is presently known that these objects have, for a 
given mass and chemical composition for the normal matter 
layers, a much lower radius than a standard WD, and also that 
they evolve in a very similar way compared to standard WDs 
(Benvenuto & Althaus 1996a, b). However, at present, it is dif- 
ficult to account for the formation of a strange quark matter 
core inside a WD star. 

In view of the above considerations, we present in this pa- 
per a detailed set of mass - radius relations for WD models 
with different assumed internal compositions, with the empha- 
sis placed on models with iron - rich composition. Despite 
the fact that many researchers have addressed the problem of 
theoretical mass - radius relations for WD of helium (Vennes 
et al. 1995; Benvenuto & Althaus 1998; Hansen & Phinney 
1998 and Driebe et al. 1998), carbon and oxygen (Koester & 
Schonberner 1986; Wood 1995 amongst others), we judge it to 
be worthwhile to extend our computations to the case of models 
with these compositions in the interests of presenting an homo- 
geneous sequence of mass - radius relations. In particular, we 
shall consider the internal layers as made up by helium (*He), 
carbon (^^C), oxygen C^O), silicon (^^Si) and iron (^^Fe), sur- 
rounded by a helium layer with a thickness of 10^^ Af* (where 
A/, is the stellar mass). We considered models with an outer- 
most hydrogen layer of 10~^ M^, (3 x 10^"* A/* in the case 
of helium core models) and also models without any hydrogen 
envelope. In doing so, we shall employ a full stellar evolution 
code, updated in order to compute the properties of iron - rich, 
degenerate plasmas properly. 

This paper is organized as follows. In Sect. ^, we present 
the general structure of the computer code that we have em- 
ployed and the main improvements we have incorporated in it. 
In Sect. ^ we describe the strategy employed in the computa- 



tions, and the numerical results. Finally, in Sect. ^ we discuss 
the main implications of our results. 

2. The computer code 

The WD evolutionary code we employed in this study is fully 
described in Althaus & Benvenuto (1997, 1998), and we re- 
fer the reader to those works for a general description. Briefly, 
the code is based on the technique developed by Kippenhahn 
et al. (1967) for calculating stellar evolution, and it includes a 
detailed and updated constitutive physics appropriate to WD 
stars. In particular, the EOS for the low - density regime is 
that of Saumon et al. (1995) for hydrogen and helium plasmas, 
whilst the treatment for the completely ionized, high - density 
regime includes ionic contributions, coulomb interactions, par- 
tially degenerate electrons, electron exchange and Thomas - 
Fermi contributions at finite temperature. Radiative opacitites 
for the high - temperature regime (T > 6000 K) with metal- 
licity Z = are those of OPAL (Iglesias & Rogers 1993), 
whilst for lower temperatures we use the Alexander & Fergu- 
son (1994) molecular opacities. 

High - density conductive opacities and the various mech- 
anisms of neutrinos emission for different chemical composi- 
tion (4He, 12c, 160, 20Ne, 24Mg, 28Si, 32s, 40ca and 56pg) 

are taken from the works of Itoh and collaborators (see Al- 
thaus & Benvenuto 1997 for details). In addition to this, we 
include conductive opacities and Bremsstrahlung neutrinos for 
the crystalline lattice phase following Itoh et al. (1984a) and 
Itoh et al. (1984b; see also erratum Itoh et al. 1987), respec- 
tively. The latter becomes relevant for WD models with iron 
core since these models begin to develop a crystalline core at 
high luminosities (up to two orders of magnitude higher than 
the luminosity at which a carbon - oxygen WD of the same 
mass begins to crystallize). With respect to the energy transport 
by convection, for the sake of simplicity, we adopt the mixing 
length prescription usually employed in most WD studies. This 
choice has no effect on the radius of the models. Finally, we 
consider the release of latent heat during crystallization in the 
same way as in Benvenuto & Althaus (1997). 

As in our previous works on WDs, we started the com- 
putations from initial models at a far higher luminosity than 
that corresponding to the most luminous models considered as 
meaningful in this paper The procedure we follow to construct 
the initial models of different stellar masses and internal chemi- 
cal composition is based on an artificial evolutionary procedure 
described in our previous papers cited above. In particular, to 
produce luminous enough initial models, we considered an ar- 
tificial energy release. After such "heating", models experience 
a transitory relaxation to the desired WD structure. Obviously, 
the initial evolution of our WD models is affected by this proce- 
dure but, for the range of luminosity and Toff values considered 
in this paper this is no longer relevant (see below) and our mass 
- radius relations are completely meaningful. 
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3. Numerical results 

In order to compute accurate mass - radius relations, we 
evolved WD models with masses ranging from 0.15 Mq to 0.5 
Mq at intervals of 5% for helium core WDs; from 0.45 Mq to 
1.2 Mq at intervals of 0.01 Mq for carbon, oxygen and silicon 
cores; and finally from 0.45 Mq to 1.0 Mq for the case of an 
iron core at intervals of 0.01 Mq. The evolutionary sequences 
were computed down to log L/Lq= -5. Mass - radius relations 
for interior compositions of ^^C, ^^O, ^*Si and ^^Fe are pre- 
sented for Tcff values ranging from Toff= 5000K to 55000 K 
with steps of 10000 K and from 70000 to 145000 K with steps 
of 15000 K. For the case of heUum WD models we considered 
Tcff values from rcff= 4000 K to 20000 K with steps of 4000 K. 
To explore the sensitivity of our results to a hydrogen envelope, 
we considered two values: Mh/M* = 10~^ (3 x 10~* M* in 
the case of He core models) and Mh/M* = 0. For the sake of 
comparison, for each of the considered core compositions we 
have also computed the zero - temperature HS models. 

It is worth noting that all the models included in the present 
work have densities below the neutronization threshold for 
each chemical composition (1.37 x 10^^ g cm~^ for helium, 
3.90 X 10^° g cm"^ for carbon, 1.90 x 10^" g cm~^ for oxy- 
gen, 1.97 X 10^ g cm""^ for silicon, and 1.14 x 10^ g cm""^ for 
iron). Such densities represent the end of the WD sequences be- 
cause electron capture softens the EOS and the stellar structure 
becomes unstable against gravitational collapse (see Shapiro & 
Teukolsky 1983 for further details). 

In recent years, both observational (Marsh 1995; Moran et 
al. 1997; Landsman et al. 1997; Edmonds et al. 1999 amongst 
others) and theoretical (Althaus & Benvenuto 1997, Benvenuto 
& Althaus 1998, Hansen & Phinney 1998, Driebe et al. 1998) 
efforts have been devoted to the study of hehum WDs. It is now 
accepted that these objects would be the result of the evolu- 
tion of certain binary systems, in which mass transfer episodes 
would lead to the formation of hehum degenerates within a 
Hubble time (see, e.g., Iben & Tutukov 1986; Alberts et al 
1996; Ergma & Sarna 1996). Connected with the age deter- 
mination for millisecond pulsars from WD cooling is the exis- 
tence or not of hydrogen flashes in helium WDs. In particular, 
detailed calculations predict that hydrogen flashes do not oc- 
cur on WDs of mass less than 0.2 Mq (see also Driebe et al. 
1998), but instead such low - mass helium degenerates experi- 
ence long - lasting phases of hydrogen burning (but see Sarna 
et al. 1998). 

With regard to the main topic of the present work, it is 
worth mentioning that Vennes, Fontaine & Brassard (1995) 
presented a set of static mass - radius relations for hot WDs 
However, the authors considered a linear relation between the 
internal luminosity and the mass, thus avoiding the computa- 
tion of evolutionary sequences. This approximation is equiva- 
lent to neglecting neutrino emission, which is not a good as- 
sumption for their hottest models. 

In a recent paper, Driebe et al. (1998) have computed the 
evolution of low mass stars from the main sequence up to the 
stage of helium WD. In that work the binary evolution has been 



mimicked by applying, at appropriate positions, large mass 
loss rates from a single star. More importantly, diffusion was 
neglected throughout the entire evolution. In this connection, 
gravitationally induced diffusion is expected to lead to notice- 
able changes in the surface gravity of their helium WD models, 
the envelope of which at the end of mass loss phase is a mixture 
of heUum and hydrogen. Indeed, during their evolution, WDs 
should modify the outer layers chemical composition making 
essentially the bulk of the hydrogen float to the surface and the 
helium sink out of surface layers. In this way, this effect causes 
the outer layers composition to approach to pure composition 
layers, the case we assumed in the present paper. Prehminary 
results to be presented below indicate this to be the actual case, 
as we suggested previously (Benvenuto & Althaus 1999). 
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Fig. 1. Surface gravity versus Tefffor 0.3 (upper curves) and 
0. 195 (lower curves) Mq helium WD models. Dotted lines cor- 
respond to models calculated by Driebe et al. (1998). Dashed 
lines represent our models with an envelope mass of 1.2 
xIO^'^Mq (with a hydrogen content by mass of Xh= 0.538) 
and 6 xIQ-^Mq (with Xh= 0.7) for flie 0.195 and 0.3 Mq 
models, respectively. Solid line corresponds to the case when 
diffusion is included in our 0.195 Mq model with an envelope 
of 1.2 X IQ-^Mq and an initial Xh= 0.538, while dot - dashed 
lines and dot - dot - dashed lines are for our 0.195 Mq mod- 
els with pure hydrogen envelopes with mass of 6 x lO"** and 
1.2 X 10~^ Mq, respectively. The location of the WD com- 
panion to the milUsecond pulsar PSR J1012 + 5307 according 
to van Kerkwijk et al. (1996) and Caflanan et al. (1998) deter- 
minations (upper and lower square, respectively) are also indi- 
cated. 

To address the problem of diffusion in helium WDs, we 
have developed a code which solves the equations describ- 
ing gravitational settling and chemical and thermal diffusion. 
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Here, we present some details of our code, deferring a thor- 
ough description to a further publication. In broad outline, we 
have solved the diffusion and heat flow equations presented by 
Burgers (1969) for the case of a multicomponent medium ap- 
propriate for the case we are studying here (see also Much- 
more 1984 for an application of the set of Burgers's equations 
to the study of diffusion in WDs). The resistence coefficients 
are from Paquette et al. (1986). To solve the continuity equa- 
tion we have generaUzed the semi - impUcit finite difference 
method presented by Iben & McDonald (1985) to include the 
effects of thermal diffusion. We have followed the evolution of 
file isotopes ^H, ^He, ^He, ^^c and ^^^O. The diffusion code has 
been coupled to our evolutionary code to foUow the chemical 
evolution of our models self - consistently. 
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Fig. 2. Evolution of the hydrogen (soUd Unes) and heUum (dot- 
ted lines) profiles as a function of the outer mass fraction q for 
the 0.195 Mq hehum WD model. Starting out from a model 
with an initially homogenous envelope with Xh= 0.538 (curves 
a), following models (b, c, d, e and f) correspond to evolution- 
ary stages characterized by log T^s= 4.096, 4.063, 4.0, 3.932 
and 3.762. 

Let us first compare our models with those of Driebe et 
al. (1998) in the case when diffusion is neglected. In Fig. 1 
we show the surface gravity in terms of Tefffor 0.195 and 0.3 
Mq helium WD models. In order to make a direct compari- 
son with Driebe et al.'s predictions, we have adopted for these 
models the same envelope mass and hydrogen surface abun- 
dance as quoted by these authors. The initial models were gen- 
erated in the same fashion as described previously. Despite the 
assertions by Driebe et al., note that our gravity values after 
the relaxation phase of our models are very similar to those 
predicted by these authors. We should remark that their "con- 
tracting models" are very different from our initial ones. In fact. 



they start with a homogeneous main sequence model in which 
nuclear energy release has been suppressed. Then, it is not sur- 
prising that they get contracting models with gravities compa- 
rable to those obtained with evolutionary models only when 
they are very cool (at TeffW 3000K for a 0.2 Mq model). On 
the contrary, in our previous works on helium WDs with hy- 
drogen envelopes, we generated our initial models from a cool 
hehum WD model, adding to it an artificial energy release up 
to the moment in which the model is very luminous. Then, we 
switch it off smoothly, getting a model very close to the cool- 
ing branch. Thus, notwithstanding Driebe et al. comments, our 
artificial procedure gives rise to mass - radius relations in good 
agreement with those found with a fully evolutionary compu- 
tation of the stages previous to the WD phase. A further com- 
parison performed with low - mass helium WD models cal- 
culated by Hansen & Phinney (1998) with thick hydrogen en- 
velopes reinforces our assertion. However, for more massive 
models some divergences appear between our results and those 
of Hansen & Phinney. Such differences are the result of the fact 
that Hansen & Phinney massive models do not converge to the 
HS predictions for zero temperature configurations, a limit to 
which our models tend. 

Now, let us consider what happens when diffusion is con- 
sidered. To this end, we have computed the evolution of a 0. 195 
Mq helium WD model with an envelope characterized by a 
mass fraction Menv/Af* = 6 x 10~^ and an abundance by mass 
of hydrogen (Xr) of 0.538, as quoted in Driebe et al. (1998). 
We begin by examining Fig. 2 in which we show the evolution 
of the hydrogen and helium profiles as a function of the inter- 
nal mass fraction for various selected values of T^a. Even in 
this case of low surface gravity, diffusion proceeds in a very 
short timescale, giving rise to pure hydrogen outermost layers. 
If we start out the computations at Tcff « 4. 1, we find that at the 
Tefi value corresponding to the WD companion to PSR J 1012 
+ 5307, our model is characterized by a pure hydrogen enve- 
lope of Mcnv/M^.. fti 4 X 10^^ (curve e). Needless to say, this 
will affect the surface gravity as compared with the case when 
diffusion is neglected. Thus, in order to accurately estimate the 
mass of that WD we do need to account for the diffusion pro- 
cess. This expectation is borne out by Fig.l, in which we have 
included the results corresponding to the situation when diffu- 
sion is included (solid line) and to the case of the model with 
assumed pure hydrogen outer layers (dot dashed lines) through- 
out its entire evolution (i.e. the conditions we would have it 
were diffusion instantaneous). In both situations we assumed 
that the total initial amount of hydrogen is the same. The dif- 
ferences in the value of the surface gravity compared with the 
case of no diffusion are noticeable. Finally, note that the track 
asymptotically merges the corresponding to complete separa- 
tion of hydrogen and helium, the structures we assumed in our 
previous works. These results clearly justifies the assumptions 
we made in our referred papers. We should also note that in 
the case of a somewhat thicker hydrogen envelope, hydrogen 
burning increases significantly, thus making the evolution con- 
siderably slower. Thus, in the plane surface gravity versus Tes, 
the asymptotic conditions of total separation of hydrogen and 
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helium would be reached far earlier in the evolution. Thus, dif- 
fusion is a fundamental ingredient if we want a solid surface 
gravity versus Tegrelation. 
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Fig. 3. The mass - radius relation for WD stars with a helium 
core. Solid and short dashed lines correspond to the cases of 
objects with an outermost pure hydrogen layer of 3 x 10~* M, 
and to models without hydrogen layer, respectively. Medium 
dashed hne represent the mass - radius relation for HS homo- 
geneous helium models. Finite temperature models are ordered 
from bottom to top with increasing Tes correspnding to (in 10^ 
K)of4, 8, 12, 16, and 20. 

In Fig. 3, we show the mass - radius relation for helium 
core models. We considered models with masses up to 0.5 Mq 
because higher mass objects should be able to ignite helium 
during previous evolutionary stages and should not end their 
lives as helium WDs. As it is well known, models have a larger 
radius the higher the Toff. It is also clearly noticeable the ef- 
fect on the stellar radius induced by the presence of an outer 
hydrogen envelope. These effects are particularly important for 
low mass models. Let us consider the case of a 0.3 Mq helium 
WD model. In the case of no hydrogen envelope, at the high- 
est Teff considered here, the object has a radius «50% larger 
than that corresponding to the HS model. If we include the hy- 
drogen envelope, the radius is «80% larger than the HS one 
(let us remind the reader that in this case we have included a 
hydrogen layer 30 times more massive than in the case of the 
other core compositions). Notice that for Teg— > 0, the radius 
of the objects tends to the HS values, as expected. In Figs. 4 
to 7 we show the results for carbon, oxygen and silicon interi- 
ors respectively. Although the effects due to finite temperature 
and the presence of an outer hydrogen envelope are also no- 
ticeable, these are not so large as in the case of the low mass 
heUum WD models shown in the previous figure. For example. 




0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 1.1 1.2 



M/M^ 

Fig. 4. The mass - radius relation for WD stars with a car- 
bon core surrounded by a helium layer with a thickness of 
10~^ M^, . Solid and short dashed Unes have the same meaning 
as in Fig. 3 but for case of the hydrogen envelope we assumed 
a mass of 10~^ M*. Medium dashed hne corresponds to the 
mass - radius relationship for homogeneous HS carbon mod- 
els. We have included in this figure the values corresponding 
to Teff (in 103 K) of 5, 15, 25, 35, 45, 55, 70, 85, 100, 115, 
130 and 145. We have also included the data for strange dwarf 
models with Ten (in lO'' K) of 10, 20, 30, 40 and 50. Notice 
that they are much more compact than standard WD models 
with the same composition and mass. In the interests of com- 
pleteness, we have extended the computations to a mass value 
of0.3MQ . 



for 1.2 Mq models, both effects are able to inflate the star only 
up to ^ 19%. This is expected because as mass increases, inter- 
nal density (and electron chemical potential ^e) also increases. 
Thus, as thermal effects enter the EOS of the degenerate gas 
as a correction oc {T / ^ef', EOS gets closer to the zero tem- 
perature behaviour, i.e. to the HS structure. As the thickness of 
the hydrogen layer is oc 51" ^ (5 is the surface gravity) it also 
tends to zero for very massive models. Note that carbon, oxy- 
gen and sihcon have a mean molecular weight per electron very 
near 2 (/z^ = 2.001299, 2.000000, 1.999364, and 1.998352 for 
hehum, carbon, oxygen and silicon respectively), thus, the dif- 
ferences in radii for a given stellar mass are almost entirely 
due to the non - ideal, corrective terms of the EOS. Because of 
this, the differences in radii are smaU, of the order of few per- 
cents. Also, in each figure we included the corresponding HS 
sequence. Notice that, for a given mass value, HS models have 
smaller radii and that there exist some minute differences even 
for the lowest T^s models. This is due simply to the presence of 
a helium (and hydrogen) layer (if present), the effect of which 
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Fig. 5. Same as Fig. 4, but for the massive carbon core models. 



on the model radius was not considered in our computation of 
the HS sequences. 
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Fig. 6. Same as Fig. 4, but for an oxygen core. 

In Fig. 4 we also included the radii corresponding to strange 
dwarfs of 0.4, 0.55 and 0.8 Mq models for T^s from reff= 
10000 K to 50000 with steps of 10000 K. These objects were 
computed assuming a carbon - oxygen composition for the nor- 
mal matter envelope, but despite the precise chemical profile, it 
is clear that they have much smaller radii than WD models of 
the same mass. 
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Fig. 7. Same as Fig. 4, but for a siUcon core. 

In Fig. 8, we show the mass - radius sequences correspond- 
ing to iron. These are noticeably different from the previously 
shown, due mainly to the higher mean molecular weight per 
electron (/Xg = 2.151344) and also to the much higher atomic 
number (Z = 26) that indicates a much strongly interacting, 
degenerate gas compared to the case of a standard composi- 
tion. 

In the case of an iron core, for a fixed mass value, the mean 
density is almost twice the corresponding to carbon and oxy- 
gen cores. Thus, it is not surprising that, for the range of Teg's 
considered here, thermal effects are less important than in the 
standard case. For example, for the 0.45 Mq iron model at 
Teff f» 25000K, thermal effects inflate its radius only by « 17%. 
For the iron core case, we have only considered models up to 
a mass value of 1.0 Mq. Higher mass objects are very near 
the mass limit for such composition (i.e. the central density be- 
comes very near the neutronization threshold, see also Koester 
& Chanmugam 1990) and would have internal densities so high 
that our description of the EOS would not be accurate enough 
for our purposes. The evolutionary sequences corresponding to 
an iron core composition are the most detailed and accurate 
computed to date, and a thorough discussion of them will be 
deferred to a separate pubhcation. 

Finally, in Fig. 9 we compare the results of our calculations 
for carbon WD models with a hydrogen envelope against the 
computations performed by Wood (1995). As we mentioned, 
the mass - radius relation have been the subject of many au- 
thors, amongst others, Koester 1978, Iben & Tutukov (1984), 
MazzitelU & D'Antona (1986), Wood (1995). Here, we shall 
compare with Wood's models since they have been thorougly 
employed by the WD community. Note that the general trend 
of our results and those of Wood is very similar. As Wood con- 
sidered more massive hydrogen envelopes (Mh/M* = 10~^) 
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Fig. 8. Same as Fig. 4, but for an iron core. We have also in- 
cluded the data corresponding to 40 Eri B, EG 50, Procyon 
B and GD 140 taken from Provencal, et al. (1998). Note the 
change in the vertical scale compared to the previous figures. 
For further discussion, see text. 



than we did, we have recomputed models with 0.6, 0.7 and 0.8 
Mq and with that hydrogen mass. We find that, for the same 
value of Mh / , Wood models have radii a bit larger than ours 
(of course, for the case of Mh/M* = 10~^ the differences are 
much larger). Note that as models evolve, such differences be- 
come smaller. 



4. Discussion and conclusions 

In this work, we have computed accurate and detailed mass - 
radius relations for white dwarf (WD) stars with different core 
chemical compositions. In particular we have considered inte- 
riors made up by helium, carbon, oxygen, silicon and iron sur- 
rounded by a helium layer containing 1 % of the stellar mass. 
With regard to the presence of a hydrogen envelope, we have 
considered two extreme values: Mh/M» = 10~^ (Mh/M* = 
3 X 10-"^ for helium core models) and Mh/M* = 0. The 
first three interior compositions are standard (according to stel- 
lar evolution theory), whereas iron - rich interiors has recently 
been suggested on the basis of new parallax determinations for 
some objects (Provencal et al. 1998). 

For computing each sequence we employed a full stellar 
evolutionary code which incorporates most of the currently 
physical processes considered relevant to the physics of WDs. 
We computed a set of evolutionary sequences for each consid- 
ered core composition by employing a small step in the stellar 
mass. We believe that these calculations may be valuable for 
the interpretation of future observations of this type of WDs. 



0.60 0.65 0.70 0.75 0.80 

M/Mq 

Fig. 9. A comparison of our results with those of Wood (1995). 
Solid lines indicate the previously presented results for carbon 
core models including an outer hydrogen envelope, solid dots 
indicate the results of Wood. As Wood considers more massive 
hydrogen layers (Mh/M^ = 10^"') than the one employed 
in most of this work, we have computed also models with the 
same hydrogen mass envelope. The results corresponding to 
these models are denoted with open diamonds. 

We have also investigated the effect of gravitational settling 
and chemical and thermal diffusion on low - mass helium WDs 
with envelopes made up of a mixture of hydrogen and helium. 
To this end, we included in our evolutionary code a set of rou- 
tines which solve the diffusion and heat flow equations for a 
multicomponent medium. For the case analysed in this paper, 
we found that diffusion gives rise to appreciable changes in the 
theoretical mass - radius relation, as compared with the case 
when diffusion is not considered (Driebe et al. 1998). 

In Figs. 4 - 8 we included the data for 40 Eri B (rcft=16700 
K), EG 50 (rcff=21000 K), Procyon B (Toff =8688 K) and GD 
140 (reff=21700 K) taken from Provencal et al. (1998). In the 
case of 40 Eri B for instance, the observed mass, radius and Tes 
are consistent with models having a carbon, oxygen and silicon 
interior and thin hydrogen envelopes. It is worth noting that 
the observed determinations are also consistent with iron core 
models with hydrogen envelope composition but for models 
with » 55000 K (models without hydrogen envelope would 
need to be even hotter). Since this temperature is far larger than 
the observed one we should discard an iron core for this object. 

On the contrary, for the cases of the other considered ob- 
jects, they fall clearly below the standard composition se- 
quences, indicating a denser interior. If we assume GD 140 
and Procyon B to have an iron core, we find that they faU on 
a sequence of a Tes compatible with the observed value. Nev- 
ertheless, the EG 50 mean radius is smaller than predicted for 
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an iron core object for the observed Teff . Thus, on the basis of 
the current observational determinations for EG 50, this WD 
seems to be even denser than an iron WD. 

At present, it seems that the physics that determines the ra- 
dius of a WD star is fairly well understood, thus the indication 
of an iron core should not be expected to be due to some error 
in the treatment of equation of state of a degenerate plasma. Ac- 
cordingly, if observations are confirmed to be accurate enough, 
we should seriously consider some physical process capable to 
produce an iron core for such low mass objects. 

Detailed tabulations of the results presented in this paper 
are available upon request from the authors at their e-mail ad- 
dress. 
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